Orientational field-dependence of low-lying excitations in mixed state of 

unconventional superconductors 
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Orientational field-dependence of the zero energy density of states (ZEDOS) is calculated for 
superconductors with the polar state (line node), axial state (point node) and 3D d-wave state. 
Depending on the gap topology and relative field direction the field dependencies of ZEDOS sen- 
sitively differ, providing us a useful and practical method to identify the gap topology. It is also 
demonstrated that for d-wave state the field rotation in the basal plane shows a sizable oscillation 
(~3%) of ZEDOS. This is directly measurable in low-T specific heat experiment in the mixed state. 

PACS numbers: 74.25.Bt, 74.25. Op 
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Much attention has been focused on various unconven- 
tional superconductors, ranging from high Tc cuprates, 
heavy Fermion materials, Boro-carbides, MgB2, to Skut- 
terudite PrOs4Sbi2. The unconventionality is associated 
with the gap anisotropy of the orbital function in addi- 
tion to the spin structure of the Cooper pair. It is quite 
important to determine the detailed nodal topology of 
the gap function; either point or line node and their lo- 
cation on the Fermi surface. The determination of these 
characteristics is expected to lead to an understanding to 
the pairing mechanism of exotic superconductors Q . 

There are several experimental methods to probe the 
gap anisotropy. One can basically distinguish the line and 
point nodes because these give rise to different and dis- 
tinct power law temperature (T) dependence in various 
physical quantities. As for the orientation of these nodes, 
there is only a few ways to probe it. The field-dependent 
thermal conductivity k{H) and the polarization depen- 
dent sound attenuation are typical ones. By measuring 
K{H,a) for different field direction a one can detect the 
location of the node in principle because the nodal quasi- 
particles (QP) under H with zero-energy transport heat 
current. In fact, a series of experiments by Izawa et al. 
have determined the location of nodes in several systems. 
These transport measurements are, however, inevitably 
involved by the scattering time effect and localization ef- 
fect of nodal QP, which hamper the determination of the 
nodal direction in some cases 0- 

Here we propose another method based on thermody- 
namics: The Sommerfeld coefficient 7 of the T-linear spe- 
cific heat at lower T is most fundamental physical quan- 
tity in Fermionic systems of interest. Since the nodal 
QP created around a vortex core in the mixed state 
sensitively reflect its gap structure, the angle-dependent 
7(ff, a) can yield characteristic oscillation pattern rela- 
tive to the nodal position under a fixed H . Recent angle- 
resolved j{H, a) measurement on YNi2B2C by Park et al 
Q demonstrates a fourfold oscillation in the basal plane 
whose amplitude 4%, nicely coinciding with k{H, a) 
experiment by Izawa etal. 0. They agree with the nodal 
direction ([100]), but disagree with the topology (point 



(line) in the latter (former)). The observed oscillation 
amplitude (~ 4%) in ^{H, a) is far off the theoretical 
prediction (~ 30%) based on the so-called Doppler shift 
argument [1^ which is qualitative in nature (see discus- 
sion in Ref. [3). In this paper we calculate the zero- 
energy density of states (ZEDOS) in the mixed state for 
various situations; its direction-dependence for the polar 
and axial state and the angular dependence for d-wave 
state, in order to examine the experimental feasibility. 
We develop a full three dimensional (3D) computation 
based on quasi-classical framework which is valid for su- 
perconductors with kpS, 3> 1 (fci? Fermi wave number 
and ^ the coherence length). This kind of calculations 
gives quantitatively reliable results. 

Anisotropic pairing is routinely analyzed within 
the separable model of pairing potential V{k, k') — 
Vbr2(fc)r2(fc')- Then order parameter takes the following 
form: A(r,fc) = '^{r)Vl{k). In the clean limit quasiclas- 
sical equations read as 
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Here hujn — 7rT(2n -I- 1) with integer n are Matsubara 
frequencies, v is Fermi velocity, $0 is fiux quantum, and 
/, f^, g are Green's functions integrated over energy 
normalized so that //^ + .9^ = 1- Fermi surface is as- 
sumed to be sphere. Order parameter '^(r) and vector- 
potential A{r) are obtained selfconsistently from the fol- 
lowing equations 
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Average over Fermi surface is denoted as (...). First, 
polar state with a line node 0) = -\/3cos0 and axial 
state with point nodes f2((/', 9) — ^3/2 sin are analyzed. 
Polar and azimuthal angle refer to the coordinate system 
with z-axis that coincide with c crystal direction. Fac- 
tors \/3 and ^3/2 assure that average of |ri((5i,0)p over 
spherical Fermi surface is unity. We are interested in 
Green's function g{r,v,uj) that describe QP excitations 
associated with vortices. The QP density of states N{E) 
with energy E relative to Fermi level is defined as 

where iVo is ZEDOS in the normal state. Green's func- 
tion g{r) spatially averaged over vortex lattice unit cell 
is denoted as g{r). We focus on ZEDOS at low temper- 
atures. This is because low temperature specific heat is 
Cs = l{B)T = 2n^h^N{E = 0, B)r/3. Therefore the 
equations are solved for T = Q.lTc- It is sufficient to 
know the Green's function only in the vortex lattice unit 
cell which is divided in mesh 41 x 41. Once the order 
parameter and vector potential are obtained selfconsis- 
tently Eqs. and (01 are solved again for u; — > 0+. 
Typically we choose uj — O.OOlirTc/h. Method of solu- 
tion is extended from and details will be described 
elsewhere 0. Here we present the results. 
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FIG. 1: a) Field dependence of ZEDOS for polar state. Scal- 
ing factor Hc2 is different for each direction. The best fit 
to low field dependence is N{E = 0,-B)/iVo = ^{B)/^n ~ 
l.U-{B/H^2)°-^^ for H II c (full circles) and 1.06- (B///c2)° *^ 
for H A- c (empty circles) . b) ZEDOS against induction B (in 
dimensionless units). 

In Fig. n field dependence of ZEDOS for polar state 
with a line node is shown ioi H \\ c (full circles) and 
H ±c (empty circles). Ratio N{E = 0, B)/No at low T 
is equal to ■^{B)/^n, where "i{B)T (jnT) is low T spe- 
cific heat of superconducting (normal) phase. We discuss 
the power-law exponent of B dependence of ZEDOS. It 
is difficult to fit the data with a single power-law func- 
tion {B/Hc2)^- At least in low field we can estimate 
N{E = 0,B)/No - {B/Ha2)°-^^ for H \\ c, i.e. very 
steep increase with field. Here is the explanation. The 



most important contribution to ZEDOS is coming from 
QP that fiow in the plane perpendicular to the applied 
field. For H \\ c geometry, those QP experience zero en- 
ergy gap. They are easily excited and extended outside of 
the vortex core even in low field in comparison to s-wave 
superconductors. The outcome is steep increase of ZE- 
DOS with field. Experimentally similar small exponent 
is observed in MgB2. The physics is analogous to the 
case of polar state. Small exponent is coming from the 
small gap at the 7r-band 8] in MgB2 while coming from 
the line node in the polar state. For perpendicular orien- 
tation H J- c (empty circles) the problem is analogous to 
that of two dimensional (2D) d-wave case for fields along 
c-axis. Power law with exponent /3 ~ 0.45 calculated 
here should be compared with self-consistent calculation 
on cylindrical Fermi surface and 2D d-wave gap function 
which reveals N{E = 0, B)/No ^ (B/i?c2)° '*^ power law 
,:j6]. For this geometry H _L c, QP in plane _L H experi- 
ence zero gap only if their momentum is in basal plane. 
Therefore, they are more difficult to excite compared to 
parallel geometry H \\ hence the exponent (3 is bigger. 

It is important to emphasize that ZEDOS in Fig. QJl) 
is plotted against B/Hc2, where Hc2 is different for each 
geometry. For polar state there is a large anisotropy of 
upper critical field H^'j ~ 2i/^. When plotted versus 
magnetic field. Fig. QJd), ZEDOS lines crosses at some 
critical field Bcr- Therefore by rotating magnetic field 
from c-axis toward the basal plane ZEDOS may increase 
or decrease depending on field value. 
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FIG. 2: a) Field dependence of ZEDOS for axial state. Scaling 
factor Hc2 is different for each direction. The best fit to low 
field dependence is N{E = 0,B)/No = 7(B)/7jv ~ 1.25 ■ 
{B/Hc2)°-^'^ for H ± c (empty circles), b) ZEDOS against 
induction B (in dimensionless units). 

For axial state with two point nodes and H \\ c, en- 
ergy gap is small only for small fraction of QP that fiow 
along the c-axis, which makes small contribution to to- 
tal ZEDOS. In this geometry ZEDOS resembles that in 
s-wave superconductors. Most of the low energy QP 
are trapped at vortex cores, at least in low field, thus 
N{E = 0,B)/Nq - B/Hc2- This is confirmed by nu- 
merical calculation shown in Fig. |2t) (full circles). For 
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FIG. 3: Angle resolved ZEDOS averaged over angle 0' in 
plane perpendicular to field, {N{e')} /No = (l/27r) / N{E = 
0, v) d(f>' for antinode (full circles) and node (empty circles) 
field direction. Line node is schematically presented with 
dashed line in the inset. Magnetic induction B = 0.0217-ffSf'''' 

rnode / rjantinode 
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field H J- c (empty circles) power-law exponent is smaller 
than that in line node polar state for both field geome- 
tries, Fig^). Roughly speaking, the larger is the angu- 
lar area of suppressed gap, the faster ZEDOS is increas- 
ing with induction B. Note that the field dependence 
N{E = 0, B)/No - {B/H^2) HB/Hc2) predicted by the 
Doppler shift calculation for point node case is far 
off the present result in Fig. |21 warning us its validity. 
Similar to the polar state, ZEDOS curves for two field di- 
rections cross at some critical field Bcr as shown in Fig. 

12b). 

We show the importance of the mutual arrangement 
of line node and magnetic field on ZEDOS. In this sense 
it is interesting to examine 3D version of — ky sym- 
metry of the gap function which is given by: 9) = 
.y/l5/4sin^ 9cos2(j). The form of the gap function is a 
natural choice for spherical Fermi surface. The question 
is, can we guess for which field direction in the basal 
plane ZEDOS is maximum based on the calculation for 
polar state? Namely, if the magnetic field H ± c is 
along the node, then all QP flowing perpendicular to 
the field experience zero energy gap (see inset in Fig. 
O analogous to polar state with H \\ c. If the field is 
along the antinode direction then QP that flow perpen- 
dicular to the field experience zero gap only if their mo- 
menta are along the c-axis, analogous to polar state and 
H Jl c. These simple qualitative arguments, suggesting 
N{E = 0, antinode) < N{E = 0, node), are misleading 
since our calculation gives opposite result. It is because 
one must take into account the contribution from QP 
that are flowing at some angle with respect to the field 
in this 3D problem. It is instructive to see how angle 
resolved ZEDOS N{E = 0,v), averaged over angle (p' 
in plane perpendicular to field, changes with angle 6' 
between QP velocity v and vortex axis. We plot this 
quantity, multiplied with weighting factor sin^' in Fig. |3| 
at some very low field. Total ZEDOS is area under the 




FIG. 4: a) Field dependence of ZEDOS for antinode (full cir- 
cles) and node (empty circles) field direction. The best fit to 
low field dependence is N{E = 0, B)/No ~ 1.14- (S/Hc2)°'^^'' 



(antinode) and N{E = 0, B)/No ~ 1.07- (B/H^; 



(node). 



b) ZEDOS is plotted versus induction B (in dimensionless 
units). In the inset ratio R = N (antinode) / N {node) is shown 
as a function of induction; denotes upper critical field 
along node direction. 



curve. Two field directions perpendicular to the c-axis 
are considered, node and antinode. For QP that flow 
perpendicular to the magnetic field, 9' — 90°, ZEDOS 
for node field direction is about 3 times larger than for 
antinode field direction. However for antinode direction 
QP experience gap node as long as the angle between QP 
velocity and field direction is 45° < 9' < 90°. In spite 
of the sin 9' weighting factor there is a significant contri- 
bution to the total ZEDOS that are coming from these 
QP. As a result at low field N{E — 0, antinode) is larger 
than N{E — 0, node) by a few percents. 

In Fig. 0]i) field dependence of ZEDOS for node and 
antinode field direction is shown. Exponent /3 w 0.32 
for antinode direction differs from /3 w 0.45 of polar 
state {H _L c) in Fig. ^ This difference comes from 
the different power expansion of the gap function in the 
node vicinity. In the antinode case, in plane perpendic- 
ular to the field, gap function can be approximated as 
|il((/), 6')| « ilo9^ in the vicinity of gap node 9 = 0. On 
the other hand, \ fl{(j),9)\ ^ \9±tt/2\ near node 9 = ±n/2 
for polar state and H ± c. The latter case is analogous 
to 2D d-wave function and H \\ c gi ving /? w 0.45. It 
was shown by Barash and Svidzinsky [iflj that tempera- 
ture dependence of specific heat is closely related to the 
exponent n of gap function power expansion near node. 
The larger is exponent n the faster specific heat Cs is 
increasing with T/Tc. Similar qualitative arguments can 
be applied to field dependence of specific heat. 

At _B = Bcr two ZEDOS curves cross, the same as 
for polar and axial state, see Fig. ^). In the inset 
of Fig. Eb), ratio R = N [antinode) / N [node) is plot- 
ted against induction B. Physics of the crossing is very 
simple and will be explained on the 3D d-wave case. 
In the present model Fermi surface is assumed to be 
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FIG. 5: ZEDOS at B = 0.0217ii'"2° as a function of angle 
a between the applied field and antinode direction. 



sphere and upper critical field anisotropy is determined 
by the gap function. For our simple 3D d-wave case 
^node ^ ^anunode Therefore, for a fixed high field 
N{E = 0, node) > N{E — 0, antinode) because along 
the node field direction the superconductor is closer to 
the normal state and ZEDOS is closer to the normal state 
value iVo. On the other hand, for B ^ Hc2 the QP exci- 
tations probe the gap structure since the biggest contri- 
bution to the ZEDOS is coming from the delocalized QP. 
It was calculated N{E = 0, node) < N{E — 0, antinode). 

The value of crossing field Bcr depends on the Fermi 
surface model. Upper critical field is also affected by 
the Fermi surface anisotropy, and can reverse the sign 
of four- fold Hc2 oscillations in the basal plane. For ex- 
ample, in YNi2B2C the gap node is along [100], [010] 
directions, implying that those are also the directions 
of -ffc2 minima. But in borocarbides Fermi surface is 
highly anisotropic. If we accept that LuNi2B2C has simi- 
lar electronic structure as Y-borocarbide, then minimum 
of upper critical field along [110] direction [ill] implies 
the decisive role of Fermi surface on Hc2 anisotropy. 
Thus, in this case the high field inequality should be 
N{E — Q,node) < N{E = 0, antinode). Since we expect 
that Fermi surface anisotropy has no role in low-field ZE- 
DOS, then the sign of four-fold ZEDOS oscillation should 
be the same for all B, i.e. there is no crossing of two ZE- 
DOS lines. 

In Fig. ISjlow field angular dependence of ZEDOS for 
field rotating in the basal plane is shown. Four-fold oscil- 
lations is what one expects from the symmetry of the gap 
function. Angular variation is ~ 3% at low fields, which 
is measurable with present experimental techniques Q. 
Note that in the 2D d-wave case Doppler-shift calcula- 
tion ^3] estimates angular variation as large as 30%. 
Parabolic-like minimum in angular dependence of ZE- 
DOS is in contrast with cusp-like minimum in 2D d- 
wave case (and cylindrical Fermi surface) pL^ lisl H^ . 



Cusp-like minimum in thermal conductivity angular de- 
pendence is predicted in the s+g model (point node) and 
observed experimentally j2]. While shape of ZEDOS min- 
imum can rule out some forms of the gap function it can 
not provide the unique answer on the question of node 
topology (line or point). It is necessary to tilt the field 
out of the basal plane and study ZEDOS to gain addi- 
tional information. This was done by measuring thermal 
conductivity and specific heat 0) with fields rotating 
around c-axis in YNi2B2C. 

We have studied the orientational field-dependence of 
the nodal QP with zero-energy in the mixed state for 
the three representative gap functions, namely the ax- 
ial, polar and dx2^y2 states. Our computation is based 
on quasi-classical approach for 3D Fermi sphere. We 
have demonstrated that the orientational dependent and 
angle-resolved specific heat measurements are an ideal 
tool to distinguish line and point nodes and to locate the 
nodal direction free from scattering time or localization 
effects associated with transport experiments and also 
that this can be feasible in the present-day technical lim- 
itations. When conducting field-rotation experiment, it 
is important to keep the field low {B < Be,-) to probe the 
intrinsic gap structure. 
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